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Abstract - The non-Markovian master equations are derived to study quantum discord dynamics 
of two qubits coupled to a common reservoir and two independent reservoirs, respectively. We 
compare the dynamics under different parameters, such as reservoir spectra and resonant param- 
eters, at high temperature and at zero temperature. The results indicate that the dynamics at 
these two extreme temperatures share similar characters, as well as differences. 



^ ' All quantum systems interact with their surrounding 
^environments. The interactions lead to dissipation and 
I— decoherence due to a flow of information between the sys- 

tern and the environment [TJ[5] . Fortunately, reservoir en- 

^ gineering techniques can suppress decoherence in labora- 
I tories and thereby make the system avoid suffering from 
^-H the environment. Engineered reservoirs arise from many 
On physical situations, including photonic crystals [5], con- 
troUable Ohmic-like environments [J, as well as optical 
1/^ and microwave cavities [5] . Recent reservoir engineering 
techniques |3H6| aim to alter the dynamics of dissipation 
^ I and decoherence in an open quantum system by modify- 
. . ing characteristics of the environment, such as reservoir 
^ temperatures, reservoir spectra and resonant parameters. 
^> Understanding which type of environment leads to faster 
5^ or slower decoherence dynamics is essential in the choice 
of the physical system for implementing realistic quantum 
devices such as a quantum computer. Therefore, it is of 
fundamental and practical importance to study the deco- 
herence dynamics of a system in structured reservoirs. 

Decoherence can be viewed as the loss of nonclassical 
correlation of a system. Entanglement and quantum dis- 
cord are two widely used measures of nonclassical correla- 
tion. Entanglement receives much attention since it plays 
a crucial role in quantum information processing. How- 
ever, as it is discoved that the unentangled states can also 
have nonclassical [3 correlation and that the use of such 
states can improve performance in some computational 
tasks IS] , quantum discord is proposed [3] and it is defined 



E-mail: stsj sjamail . sysu.edu.cn 



as distinction between quantum and classical aspects of 
correlation in a composite quantum state. Entanglement 
and quantum discord can behave very differently under 
certain reservoirs. For example, under Markovian envi- 
ronments |10j quantum discord is more robust than the 
entanglement against decoherence; under independent jllj 
or common |12| zero-temperature non-Markovian reser- 
voirs, quantum discord presents sudden changes |13| while 
entanglement shows oscillations with or without suddden 
death. Experimentally, the above distinguished behaviors 
have been obtained [14] . Recent results suggest that a 
change in quantum discord could be used as an indica- 
tor of failure of local operations and classical communi- 
cations }15| . Besides, the quantum thermodynamics effi- 
ciency of a photo-Carnot engine in terms of quantum dis- 
cord of an atomic pair can exceed its classical value |16j . 
In this letter, we employ quantum discord to character- 
ize the quantum correlation present in a two-qubit sys- 
tem. We consider these two qubits subject to two dif- 
ferent situations, the independent and common environ- 
ments, and aim to find out under which situation coher- 
ence can be preserved longer in high-temperature and in 
low-temperature regions, respectively. With Einstein con- 
vention sum being adopted, the Hamiltonian for the case 
of independent reservoirs reads {h = 1) 



and the Hamiltonian for the case of a common reservoir 
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takes the form 

= ujao\a\ + ^c^fc4afe + {g\B + o\B^), (2) 

where uja is the transition frequency of the qubits, a\_ are 
the system raising and lowering operators of the zth qubit, 
and J5' = Y.k 9k°'\ = Efc 9kak) with gl (gk) being the 
couphng constants. In the hmit of a continuum of reser- 
voir modes l^fep — )• / dijjj[uj), where J(a;) is the spec- 
tral density function, charactering the reservoir spectrum. 
Here the index k labels the reservoir field modes with fre- 
quencies (wfe), and a\ and a\. {a\ and ak) are their 
usual creation and annihilation operators, respectively. 

Let us now proceed to the master equation for this two- 
qubit system using the method suggested in [T7]- In the 
weak coupling limit, assuming an initially factorized state 
and a thermal reservoir, we obtain from eq. ((T]) a secularly 
approximated non-Markovian master equation in the in- 
teraction picture 

dp . -J 
— = - 4:Kip - lK2JoP 
at 

+ 2{ki + ^ii)K-p + 2(ki - fii)K+p - AjJiKop (3) 
and stating from eq. we have 



dp 

— = - Akip - 2ki Ji/9 - 2ifi2J2P 
at 

+ 2(ki + pi)J-p + 2(ki - pi)J+P - in2JoP 

+ 2(ki -I- t-ii)K-p + 2(ki - Hi)K+p - 4/ii/voP (4) 

where Jo,i,2,± and Kq ± are superoperators defined as 







- pcri^' - pcr^^\ 


(5) 
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J+P = cr+Vcr_' 
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, (2) (2) 


(10) 




K^P - ^^V-^i^^ 


, (2) (2) 


(11) 



Kqp 



(1) (1) , (1) (1) 

p + pa\ a_ — p 



(2) (2) , (2) (2) 



(12) 



and the time dependent coefficients Ki(t) and pi(t) (i = 
1,2) can be expressed as power series in the system reser- 
voir coupling constant a [to be defined in eq. (fT7)) ]. For 
weak coupling (i.e. when a ^ 1), one can stop the ex- 
pansion to the second order and obtain analytic solutions 



for these coefficients. In thermal equilibrium, these coeffi- 
cients read 

Ki{t) = ^J dr J dujJ{uj){l + 2N{uj))coa{uj-uja)T, (13) 



P2{t) - - 



l + 2N{u))sin{u~Ua)T, (14) 



dr / dujj{uj) cos{uj — uJa)T, (15) 



dr / dujJ{uj) sin(a; — OJa)^ 



(16) 



where N{uj) = (e'^/*-'^^ — 1)~^ is the average number distri- 
bution of the reservoir thermal excitation with frequency 
io at the initial time, ks the Boltzmann constant, and T 
the reservoir temperature. A closed form for the expres- 
sions of the time-dependent coefficients can be obtained 
in the high-temperature and zero-temperature limits, i.e., 
for l + 2N{uj) « 2kBT/u} and \ + 2N{uj) « 1, respectively. 
The solutions to the master equations ([3]) and @ can be 
obtained in the supplementary material. 

All the dynamics characters of the system are deter- 
mined by the time dependent coefficients, K2, p,2, ki and 
^1, in the master equations ([3]) and The difference 
between eq. ([3]) and eq. ^ is that the latter induces su- 
peroperators Ji,2.±- From eq. ~ cq. (jl2p . we know that 
Ji,2.± are caused by the coupling between qubits, while Jq 
and Kq^± are introduced by each independent qubit. Now 
two questions come up; do superoperators ^i,2,± speed up 
or slow down the decoherence? does the role of Ji,2,± vary 
between high and low temperature region, between reso- 
nant and off-resonant region, between different spectra of 
reservoir? We will answer these questions in the following. 

In order to compare the dynamics for different types of 
reservoirs, we consider a class of spectral densities. The 
spectral densities we examine are of the form 



J(c.) 



2 1- 
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(17) 



which is classified as sub-Ohmic if < s < 1, Ohmic if 
s — 1, and super-Ohmic if s > 1. We consider three ex- 
amples with s = 1/2, 1 and 3 corresponding to sub-Ohmic, 
Ohmic, and super-Ohmic spectral densities, respectively. 
The parameter Wc is the cutoff of high frequency of reser- 
voir and is connected to the reservoir correlation time 
Tc ~ ^c^- On the other hand, the dimensionless cou- 
pling constant a is related to the relaxation time scale 
Tr ~ Q;~^, over which the state of the system changes, 
in the Markovian limit of flat spectrum. Generally, the 
dynamics of a quantum system comprises three different 
dynamical effects occurring at three different respective 
time scales. Firstly, the dynamics of decoherence occurs 
at a time scale of the order of the relaxation time scale Tr, 
which is defined by the properties of the reservoir. Sec- 
ondly, nonsecular terms cause oscillations, occurring over 



a typical time scale Tq 



of the system. Finally, the 
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Fig. 1: (Colors online) Time evolutions of the quantum discord of the two qubits with an initial Bell state I'i/) — {\e1g2) + 
|<?ie2))/\/2- (a) ujc = lOoja for independent reservoirs, (b) ujc — oJa for independent reservoirs, (c) u>c = O.Sua for independent 
reservoirs, (d) ujc — lOtJa for a common reservoir, (e) ljc = for a common reservoir, (f) uic = O.ScJa for a common reservoir. 
The black solid line is for sub-Ohmic reservoir, blue dashed for Ohmic reservoir and pink dash-dotted for super-Ohmic reservoir. 
The insets show the same plots but with much larger values of horizontal axis ujat. We have set ksT / {huja) = 100. 



non-Markovian memory effects happen for a time shorter 
than or of the order of the reservoir correlation time scale 
Tc. In the present work, we focus on both short-time and 
long-time evolution at the secular regime characterized by 
the condition Tq <C tv [I] . In all the plots presented below, 
we set the coupling constant ~ O.OIwq, which satisfies 
the secular approximation condition Ta ^ r^. 

Now let us consider a bipartite state pab to review 
the definition of quantum discord [HIIIHI- The idea of 
quantum discord grows out of the fact that the quan- 
tum mutual information of the state pab may be de- 
termined in two nonequivalcnt ways. The first is ob- 
tained by I{pab) = S{pa) + S{pb) - S{pab), where 
S{p) = —Tr[plogp] is the von Neumann entropy of the 
state p and pa{B) = 2^?'b(A)[pab] are the partial traces 
over the two subsystems. On the other hand, one may ac- 
quire the quantum mutual information based on condition 
entropy, that is, Ja = max^Yl^y[S{pA) - S{pAB\{l\k})]^ 
where the maximum is taken over the set of position- 
operator-valued measurement (POVM) {Ofe} in partition 
B, S{pAB\{n.k}) = J2kPkSiPk) is the quantum condition 
entropy, pk = {I (^Uk)PAB{I(?)Uk)/TriI(^Uk)PABiI (?) 
Yik) is ti^® conditional density operator corresponding to 
the outcome labeled by fc, and pk = Tr{I (8) Ylk)Pab{I ® 
Yik)- Here / is the identity operator on subsystem A. Fi- 
nally, the quantum A discord is defined in terms of the 
mismatch Q{pab) = 1-{pab) ~ Ja{pab)- Similarly, one 
can also define the B discord through the entropy of con- 
ditional states of system B. In the following, we will em- 
ploy the measure of quantum discord to study under which 
situations the coherence can be maintained longer in high 



and low temperature region, respectively. 

We begin addressing the high-temperature region, 
where we choose a temperature such that ksT — lOOwa- 
Fig. [T] shows us the time evolution of the quantum dis- 
cord of the two qubits with an initial Bell state j^*) = 
{\s1g2) + \gie2))/ V^, coupled to a common reservoir or 
different individual reservoirs at high temperature, for 
three different values of resonance parameters, ujc = lOcJa, 
uJc — i^a and uJc = O.Swa- We first compare quantum 
discord evolution under a common reservoir with that un- 
der independent reservoirs. As shown in fig. [TJc) and (f), 
the quantum discord is preserved in a steady state in the 
case of a common reservoir, while in the case of indepen- 
dent reservoirs it decays exponentially with time. This 
phenomenon is due to the difference between eq. ^ and 
eq. (|4]), that is, the correlated shifted frequency induced 
by J2 and decay rate caused by Ji^± between the two 
qubits. Note that the behavior agrees with the results 
of [H] where the entanglement dynamics for two oscilla- 
tors or two qubits in the same environment was discussed. 
The model that a two-qubit system is coupled to a com- 
mon bath will attract much attention, because under this 
model one may preserve quantum coherence at high tem- 
perature, which makes it a major distinct feature with re- 
spect to the case of separate baths. Similar results can be 
obtained by comparing fig.[lla) with (d) and (b) with (e). 
We next compare the three different Ohmic-like reservoirs 
and find out which one induces the slowest decoherence. 
From fig. m we can see that the Ohmic reservoir induces 
the slowest decoherence, while for the super-Ohmic and 
sub-Ohmic reservoirs the discord decay in a very similar 
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Fig. 2: (Colors online) Time evolutions of the quantum discord of the two qubits with an initial Bell state \^) = {\e1g2) + 
l5ie2))/\/2. (a) ujc — IQoJa for independent reservoirs, (b) cjc = i^a for independent reservoirs, (c) ojc ~ O.Soja for independent 
reservoirs, (d) uJc = lQu}a for a common reservoir, (e) ujc = i^a for a common reservoir, (f) uic = O.Sa;^ for a common reservoir. 
The black solid line is for sub-Ohmic reservoir, blue dashed for Ohmic reservoir and pink dash-dotted for super-Ohmic reservoir. 
The insets show the same plots but with much smaller values of vertical axis quantum discord and horizontal axis uiat. We 
have set T = 0. 



manner, both faster than the Ohmic case, in short-time re- 
gion. Therefore, if one is able to modify the natural reser- 
voir spectrum into an Ohmic form, one would slow down 
decoherence with respect to the sub-Ohmic and super- 
Ohmic ones. This conclusion is consistent with the results 
found in [19j and [20[ . where a quantum harmonic oscil- 
lator and two cavity fields are investigated, respectively. 
However, it should be mentioned that fig. [IJa) also indi- 
cates that the super-Ohmic independent reservoirs induce 
the slowest decoherence in long-time region. Finally, the 
effect of resonance parameter will be reviewed. Changing 
this parameter corresponds to shifting the qubit frequency 
with respect to the reservoir spectrum. This allows us to 
control the effective coupling between the system and the 
environment. For p 1 {p = io^jiOo) the system is off- 
resonant with respect to the peak of the reservoir spec- 
trum. Accordingly, for p 3> 1 the system is resonant with 
respect to the peak of the reservoir spectrum. From pre- 
vious results [22], we expect to see different dynamics in 
the p ^ 1 and p'^ \ regimes. Comparing fig. [TJa) with 
(c), (d) with (f), we find that, whatever type of reservoir 
the system is under, for p 3> 1 the decoherence process 
is faster than that for p <C 1 in slrort-time region. This 
dues to the overlap between the frequency of the system 
and the reservoir spectrum in the resonant case. But in 
long-time region the decoherence process is significantly 
slower for p 3> 1 than for p <^ \ under the influence of 
independent super-Ohmic reservoirs, as shown in fig. [IJa) 
and (c). Thus in the p ^ 1 case the effective coupling of 
the system to the reservoir is obviously weaker than that 



in the p case for super-Ohmic reservoir in long-time 
region. In summary, the results show that (i) A common 
reservoir forms a steady state at high temperature, which 
is superior to the case of two separate reservoirs; (ii)For 
Ohmic reservoir the system shows the longest quantum 
coherence in each situation except for the situation of in- 
dependent reservoir in resonant region, where for super- 
Ohmic reservoir it appears to have the weakest effect on 
the system in long-time region; (iii) Under sub-Ohmic and 
Ohmic reservoirs, the coherence will be maintained longer 
in the off-resonant region, while for super-Ohmic reservoir 
we should consider it for two cases. In case of short-time 
region, quantum correlation can be preserved longer in the 
off-resonant region, however, in case of long-time region 
this correlation can be protected for a longer time. 

We now move on to the zero-temperature region, i.e. 
r = 0. From previous studies [22] on open quantum 
system interacting with zero-temperature reservoirs, we 
expect an apparently slower loss of coherence with re- 
spect to the r ^ case. Nevertheless questions still re- 
main. Does the zero-temperature common reservoir in- 
duce a steady state like the high-temperature common 
reservoir? Does Ohmic environment still show the weak- 
est dissipation-noise effect on the decoherence dynamics? 
Does the resonant parameter have the same influence as 
at high temperature? In order to answer these questions, 
in fig. [2] we plot time evolution with the same parameters 
as fig. [T] except for temperature. As shown in fig. [H the 
zero-temperature common reservoir not only doesn't in- 
duce a steady state, but also speeds up the decoherence 
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Fig. 3: (Colors online) Time evolutions of the quantum discord of the two qubits with an initial separable state l^*)! = \e1g2) 
or |>I')2 = \e1e2) for a zero-temperature common reservoir, (a) ljc — lOuJa with (b) cjc = i^a with (c) ojc = 0-3iOa 

with l^')!, (d) iUc = lOuJa with |*I')2, (e) lJc = coa with |*I')2, (f) i^c = O-Suja with |^')2. The black solid line is for sub-Ohmic 
reservoir, blue dashed for Ohmic reservoir and pink dash-dotted for super-Ohmic reservoir. The insets show very-long-time 
scale dynamics. We have set T = 0. 



rate, compared with the case of independent reservoir, on 
condition that other parameters are the same. As for the 
second question, we discuss short-time and long-time re- 
gion separately. From the insets of fig. [H we find that 
super-Ohmic environment shows the strongest decoher- 
ence, while the sub-Ohmic and Ohmic reservoirs decay in 
a very similar manner, both slower than the super-Ohmic 
case, in short-time region. By contrast, it is not easy to 
say which type of spectra induces the slowest decoherence 
in long-time region, since it dependents on resonant pa- 
rameters. For example, sub-Ohmic reservoir shows the 
weakest effect in case Wc = O.ScJq, but in case uJc = IOwq 
super-Ohmic reservoir does. Finally, fig.[2]also tells us that 
the effects of resonant parameters are similar to those at 
the high temperature, that is, the decoherence process is 
faster for p 3> 1 than p ^ 1 under the influence of Ohmic 
or sub-Ohmic reservoir, but conversely in case of super- 
Ohmic reservoir. All in all, the decoherence dynamics in 
zero-temperature and high-temperature share some inden- 
tical characters, but are different from each other as to 
whether a steady-state will be produced under a common 
reservoir. Furthermore, at zero temperature what type 
of reservoir leads to the weakest decoherence depends on 
resonant parameters and time regions, while at high tem- 
perature Ohmic reservoir always shows the weakest impact 
on the system. Lastly, the roles of resonant parameters are 
similar to those in high temperature region. 

As well as avoiding decoherence, quantum coherence 
must also be generated. In fact, it has been reported that 
two qubits coupled to a single mode Plj or a common 
heat bath |25| can be entangled by purely dissipative dy- 



namics. Here, we focus on creation of quantum discord 
of two qubits interacting with a zero-temperature com- 
mon bath and consider the behavior characteristics for 
varies of parameters. In fig. [3J we plot the dynamics of 
this system with different initial states = \e1g2) and 
|^')2 = \e1e2), respectively. A first look at the plots of [3] 
shows that quantum discord can be generated while the 
evolution varies with reservoir spectrum, resonant param- 
eter and initial state. For the initial state I^P)! = \e1g2), 
we find from fig. [3] (a) ^ (c) that the quantum discord 
increases and approaches to a definite value monotoni- 
cally. This conclusion is similar to the result found in |26) 
where dynamics of bath-induced entanglement was an- 
alyzed. And one may also sec that the time evolution 
of quantum discord is changed with the type of spectral 
densities for the same resonant parameter. For example, 
as shown in fig. [3] (a), the super-Ohmic environment in- 
duces the strongest quantum discord oscillation; next is 
the Ohmic case; while the sub-Ohmic environment causes 
a relatively weak quantum discord oscillation. In addition, 
for a fixed type of spectral density, the behavior of quan- 
tum discord varies with resonant parameter. Take super- 
Ohmic for example, it oscillates very quickly for p ^ 1 
region, but no oscillations are present for p ^ 1. On the 
other hand, for the initial state |\I')2 = |eie2), the quantum 
discord increases to a fixed maximum value, then decreases 
gradually without any oscillations, as shown in fig. [21(d) ^ 
(f). For this kind of initial state, the functions of resonant 
parameters depend on reservoir spectra. For sub-Ohmic 
or Ohmic reservoir, resonant parameters in p 3> 1 case will 
shorten the life of quantum discord, while in p <C 1 case 
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they prolong the Hfe of quantum discord. Instead, adverse 
results will been observed for super-Ohniic reservoir. 

In conclusion, the non-Markovian master equations are 
derived to study the decoherence dynamics of two qubits 
weakly coupled to a common and two separate bosonic 
thermal baths, respectively. We compare the dynamics at 
zero temperature and at high temperature under differ- 
ent parameters, such as resonant parameters and reservoir 
spectra, and we find that there are both similarties and 
differences between the decoherence processes at these two 
extreme temperatures: firsly, coherence of two qubits in- 
teracting with a common reservoir can be preserved longer 
than that with two independ reservoirs at high tempera- 
ture, while the reverse result will be seen at zero tempera- 
ture; secondly, on the one hand for the super-Ohmic reser- 
voir and sub-Ohmic reservoir the discord decay in a very 
similar manner and both are faster than Ohmic case in 
short-time region at high temperature, on the other hand 
the discord decay in a similar manner for sub-Ohmic and 
Ohmic reservoir and both are slower than supcr-Ohmic 
reservoir case in short-time region at zero temperature; 
thirdly, the difference between the decoherence processes 
in resonant and off-resonant region at high temperature is 
similar to that at zero temperature; finally, non-classical 
correlation of two qubits coupled to a common reservoir 
can be generated at zero temperature. The results in this 
paper not only have fundamental interest for the study of 
quantum discord, but also can be utilized to design quan- 
tum computation in the future. 
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The solutions to the master equations will be given in this supplementary 
material. We aim to demonstrate that the case of independent reservoirs is very 
different from the case of a common reservoir with the appearance of superop- 
erators Ji 2 ± ■ 



1 The master equations 



The non-Markovian master equation can be obtained by The non-Markovian 
master equation for two qubits coupled to two independent reservoirs reads 



at 

+ 2{ki + ^ll)K^p + 2(ki - pLi)K+p - 4fiiKop, (1) 

and the master equation for two qubits coupled to a common reservoir takes 
the form 



dp 

— = -4kip - 2ki Jip - i2p.2J2P 
dt 

+ 2{ki + pi)J-p + 2{ki - pi)J+p - iK2Jop 

+ 2(ki + pi)K^p + 2(ki - pi)K+p - 4piKop, (2) 

where Jo,i,2,± and Kq^± are superoperators defined as 



7^ (1) (1) I (2) (2) fo\ 

*The supplementary material is related to the manuscript "Quantum discord dynamics in 
structured reservoirs", which is submitted to Europhysics Letters. 
tE-mail:suzhikun(9!163.com 
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X+p = a(^ai^'+.fp.i^ (4) 

2 2 ' ^' 

J,p^a(^^p + af^p-paf)^pafK (6) 

T (1) (2) , (1) (2) , (1) (2) , (1) (2) 

Jip = CT_ (tX p + (T\'a_' p + p(j\'a\' + pa\'(T\' , (7) 

7 (1) (2) , (2) (1) (1) (2) (2) (1) 

J2P = (7j_ (J_ p + (j\ (J_ p — pa\'(j_' — p(7\ (J_ , (8j 

7 (1) (2) , (2) (1) /fiN 

J+P ^ <y+ pr^- + P<y- ■ (9) 

2 The solutions 

2.1 For the case of separate reservoirs 

The solution to master eq.([T|) can be obtained with the algebaraic approach[T] 
because the superoperators herein satisfy SU(2) Lie algebraic communication 
relations, i.e. 

[Jo,^±,o] =0, 

[i^o, K±] = ±A±, [K^,K+] = -2/io. (10) 
By directly integrating eq.(IT]), the formal solution is obtained as 

p{t) = e-reJo^ofg/o'<it(^oifo+^+if++^-if-)^(o)^ (H) 

where T is the time ordering operator, 

jo = -2* / dtK2{K2 + P'2),l^0 = -4^1, 

Jo 

I/+ = 2(ki = 2(ki +pi),r = 4 / dtKi. (12) 

Jo 

The exponential function of superoperator in cq. (jlip can be factorized in 
the form 

where fco,A;+ and A;_ satisfy the following differential equations 
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/e+ = — v-k^ + 1^0^+, 



(14) 



ko = vo- 2!/_fc+, 

k- = V- exp(fco), 
and using the following relations 



(15) 
(16) 



e-^" = (c/^f +a(^)s/.f)(c/.f + aW./.f) 
.p.(c/.i+.«s/.f)(c/,,f +af./.f), 



, J2) (2) 



2^(1) J2) _(1) (2) 



(17) 
(18) 



-p = p + j.ia'ypa'_^^ + a'^^paX') + (i-)VL^VL^V^V^VV^\ (19) 



• p • [1 + (e* - l)a«<7L^)]}[l + (e* - l)af aL')]. 
Therefore, one can obtain 



Pit) 



( Pllit) Pl2{t) Puit) Puit) \ 

P2l{t) P22{t) P2i{t) P2i{t) 

P3l{t) P32{t) P33(t) p34{t) 

\ P4l{t) P42{t) PAsit) P44(i) J 



(20) 



(21) 



with 



pii(t) = (e*=° + 2k+k- + e-'=°4fc2 )pn(0) + (fc+ + e-'=°4fc_)p22(0) 
+ (fc+ + e-'=°4fc_)p33(0) + e-'=''4p44(0), 

P22(i) = (k. + e-^"k+e_)pn{0) + (1 + e-'^''fc+fc_)p22(0) 
+ (1 + e-'=°/c+/c_)p33(0) + e-'=U+p44(0), 



P33(t) = (fc- + e-'=°fc+e)pii(0) + e-'=°A;+A;_p22(0) 
+ e-'=°fc+fc_p33(0) + e-'=°fc+p44(0). 
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P44(i) = e-'^^fc^ pii(O) + e^'=''fc_[p22(0) + pasW] + +e-'=V44(0), 



P2i(t) = (e"°/2 + e-'="/2A;+A:_)p2i(0) + e-^^^'^k+p^m, 
P3i(t) = (e'=°/2 + e-'=°/2A;+A:_)p3i(0) + e-'=°/2A;+p42(0), 
Pi2(t) = (e'=«/2 + e-'="/2fc+fc_)pi2(0) + e-'="/2fc+p34(0), 
P42(t) = e-'=«/2;._^g^(o) _^ e-'="/V42(0), 

Pi3(t) = (e'=°/' + e-'="/2A;+A;_)pi3(0) + e-'''°/^k+p2m, 
P43{t) = e-'=«/2fc-P2i(0) + e-'="/2p43(0), 
P24(t) = e-'^«/2^_pi3(0) + e-'="/2p24(0), 
P34(t) = e-'=''/2fc_p^2(0) + e-'="/2p34(0), 

Pl4(t) = Pl4(0),P23(t) =P23(0),P32W =P32(0),P4l(t) =P4l(0). (22) 

From eq. p^ . we can see that Jo,K± o can't affect the anti-diagonal elements 
of the density matrix. In the next section, we aim to find out wheather the 
anti-diagonal elements of the density matrix will be affected in the case of a 
common reservoir. 

2.2 For the case of a common reservoir 

The formal solution to eq.© is obtained as 

The superoperators satisfy the communication relations 

[Jq, '^±,1,2] = 0, 
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r (2) (1) (1) (1) , (1) (1) (1) (2) , (1) (2) (2) (2) , (2) (2) (2) (l),-, 
[Jl,J2] = 0, 

[Ji, J_]p = af)a^l~^pa^l~^ + a^^^a^^'^paf + a^^~^pafaf) + ai^V^^V^^), 
[J„ J+]p = -a^^afpa^ - af^a^^^pa^^ - a^pa^af^ - pai^V^^), 
[J2, J_]p = a^)a^^^paf + af^a^l^pa't' - a'^pa^^'a^ - a^^^pafa^^\ 

[J2,i^_]p = a(^V(^)paW paf - ai^V^^V^ - ai^a^^f , 

[J2, K+]p = -a^^^aPpa^^^ - afai'^pa^^ + <7«p^i'Vi') + afpa^a^^, 

t r ^ W (2) (2) (2) , (2) (1) (1) (1) (1) (1) (1) (2) (2) (2) (2) 

1/ (2) (1) (1) (1) (1) (1) (1) (2) , (1) (2) (2) (2) (2) (2) (2) (l)-, 

— cr)_ a_'a^' p — pa\' ay a\' + a\' ay a\' ay p — pa^_'aya_'ay), 



[J_,J+]p = 



(1) (1) (2) 



(1) (2) (2) (2) (1) (1) 




[J_,if_]p = 0, 




and 




(24) 



We then approximately] the exponential in the right-hand side of eg. ([25)1 . 
The first-order approximant is given by 



Solving eq. I\2i\ even for its first-order approximant is a non trivial task, but 
we can see from eq. ([25|) that the superoperators Ji,2,± will affect on the initial 
state of the system and further caculation shows that Ji,2,± will affect on all 
elements, including anti-diagonal elements, of the density matrix. It should be 
pointed out that we evaluate dynamics of qubits coupled to a common reservoir 
numerically in the manuscript. We are intertested in what will be induced 
by the difference between the case of independent reservoirs and the case of 
a common reservoir and whether it will vary from high-temperature region to 
zero-temperature region. 
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X'g/o dt(eiJi+e2J2+e+J++e-J-+uoKo+v+K++u-K-) 

Ji-l-e2 J2-l-e+ J++e- J-)g-it(i'o-ft^o+!^+-ff++i^--ff-) 
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